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Abstract 
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In this paper, a link adaptation and untrusted relay assignment (LAURA) framework for efficient and reliable 
1 wireless cooperative communications with physical layer security is proposed. Using sharp channel codes in different 

, transmission modes, reliability for the destination and security in the presence of untrusted relays (low probability of 

interception) are provided through rate and power allocation. Within this framework, several schemes are designed 

7— I 1 

^ ' for highly spectrally efficient link adaptation and relay selection, which involve different levels of complexity and 

o 

£ — ■ channel state information requirement. Analytical and simulation performance evaluation of the proposed LAURA 

o- 

\Q ' schemes are provided, which demonstrates the effectiveness of the presented designs. The results indicate that power 

00" adaptation at the source plays a critical role in spectral efficiency performance. Also, it is shown that relay selection 

O! 

based on the signal to noise ratio of the source to relays channels provides an interesting balance of performance and 

t-H ■ 

J!" 1 complexity within the proposed LAURA framework. 

X: 

Index Terms 

Amplify-and-forward relaying, cooperative communications, link adaptation, physical layer security, relay selec- 
tion, untrusted relay. 

I. Introduction 

Diversity and link adaptation are two enabling techniques to facilitate high performance communications over 
wireless fading channels. In cooperation for diversity, the relays are to assist a reliable data transmission from the 

A preliminary report on this work has been presented at the IEEE International Conference on Telecommunications, Cyprus, May 2011 |T|. 
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source to the destination. One major challenge hindering the practical adoption of cooperative wireless communica- 
tions is the security as the relaying nodes may in fact be able to eavesdrop on the source destination communications. 
In this paper, wireless link adaptation solutions are proposed to facilitate both reliability and physical layer security 
for cooperative communications in presence of untrusted relays. 

The cooperation of the source with a relay node may include a so-called service level trust, i.e., the relay node 
indeed performs its expected function as a relay in the network. However, this cooperation may not necessarily 
include a data level trust, i.e., the relay may not be supposed to extract (decode) useful information from the source 
destination communication. The information theoretic aspect of this problem is investigated in [2|-|4|. Specifically, 
an upper bound for the achievable secrecy rate in this setting is presented in |2|. 

Link adaptation by rate and power control could highly improve the performance of cooperative communica- 
tions over time-varying channels. In 0, the capacity of adaptive transmission over cooperative fading channel is 
considered for amplify-and-forward relaying, where three different adaptive techniques are investigated. In [6|, the 
performance of cooperative communications with relay selection and un-coded adaptive modulation is investigated. 
In 0, a cross-layer approach to optimize the spectral efficiency of the relay channel employing adaptive modulation 
and coding in conjunction with cooperative automatic repeat request is proposed. 

A practical code design approach to physical layer security based on capacity achieving low density parity check 
(LDPC) codes is introduced in [8 |. The setting is a wiretap binary erasure channel involving a source, a destination 
and an eavesdropper. In fact, constructive approaches for (imperfect) physical layer security may be set up to ensure 
a reliable source destination communication, while maintaining a high probability of error for the eavesdropper. To 
this end, design of sharp punctured LDPC codes for which the bit error rate (BER) curve fall sharply from high 
BERs to very low BERs (steep waterfall region), is considered in [9|. Yet, adaptive transmission may be exploited 
to enhance the system performance over time-varying channels both in terms of security and reliability [1 1. 

In this paper, a framework for cooperative communications through link adaptation with untrusted relay assignment 
(LAURA) is proposed. The purpose is to utilize the cooperation of arbitrary number of relays for reliable communi- 
cations, while ensuring they cannot decode useful information from their relayed signal. Different LAURA schemes 
are set up based on reasonable available CSI of the network that is already necessary for quality of service (QoS) 
provisioning over wireless fading channels. The resource degrees of freedom, i.e, transmission power, transmission 
rate, and cooperating relays, are exploited in an optimized manner to address the QoS reliability and security 



3 




Fig. 1. Network topology; 7 and £ indicate the channel SNR and the distance. 

requirements of the cooperative communication networks dealing with untrusted relays. Specifically, several power 
adaptation and relay selection strategies are proposed for LAURA that are designed for high spectral efficiency 
communication based on discrete rate adaptation with sharp channel codes and different levels of complexity and 
CSI requirements. Analytical and simulation performance evaluation of the proposed LAURA schemes are provided, 
which demonstrates the effectiveness of the presented designs. The results indicate that the LAURA design with 
source (only) power adaptation and relay selection based on (only) source relay channels CSI provides an interesting 
balance of performance and complexity within the proposed LAURA framework. 

The paper is organized as follows. Section [II] presents the system model and describes the problem under 
consideration. In Section [III1 the link adaptation and relay selection problem is introduced and its exact solution is 
presented. Section [IV] explores the scenarios where source and relays transmit with constant power. In Section [V] 
suboptimal relay selection strategies are investigated. Section [VI] presents the simulation and theoretical results for 
the proposed LAURA schemes. Section \VU\ concludes the paper. 

II. System Model and Problem Statement 

A wireless communication system with one source node (S), one destination node (D) and a set of JVr available 
relay nodes, {UR,-|, i g .Mr = {1, 2, j?Vr}}, is considered (Fig. [T}. There are Nq cooperating relays, denoted by 
the set A4c, that are selected from the set of available relays A4r. The cooperation protocol is amplify-and-forward 
(AF). The cooperating relays are assumed trusted at the service level and untrusted at the data level. Service level 
trust entails true CSI feedback, power adaptation according to source schedule and forwarding the amplified version 



of received signal without modification. Non-cooperating relays are also service level trusted, i.e., they are only 
capable of listening to the transmitted signal in the first half of transmission interval as they are expected. Since, 
all relays are data level untrusted, we have a security constraint for cooperating and non-cooperating relays, i.e., 
transmitter wishes to prevent all relay nodes from extracting useful information from their received signal. 

In the first phase of the transmission, the S node transmits signal x to the D node and relays. The received signal 
at the D node and at the UR, node are respectively 

Vs& = K&X + W i (1) 

Vsi = Kt x + w si , (2) 

where h S i, and h S i denote the Rayleigh fading coefficients between S and D nodes, and S and UR, nodes, 
respectively. The noise is denoted at the UR^ node as and at the D node as Wd- The Nq cooperating relays 
amplify the received signal and transmit it to the D node in the second phase of cooperation. The received signal 
at the D node from the UR 2 node is 

Via = Gi h lA y si + w' d , (3) 

where G,; is the UR^ node amplifier gain and ha and w' A are Rayleigh fading coefficients from URi node to the 
D node and noise at the D node, respectively. 

The S-D, S-URi and UR^-D channels, for 1 < i < Nr, are independent Rayleigh fading channels with SNRs of 
7sd, 7si and 7^, respectively. The SNRs are exponentially distributed with parameters -J^, J- and and, probability 
density functions of / s d(7sd), fsi{%i) and fid^a), respectively. With full CSI assumption, the instantaneous SNRs 
of the three channels are assumed known at the onset of each frame interval at the S node. Assuming maximum ratio 
combining (MRC) of signals received from S node and Nc cooperating relay nodes at the D node, the equivalent 
SNR of AF relaying protocol at the D node is iflOl 

y S d -+ 

ieMc 

where 



7eq = ^77sd + ^2 li, (4) 



f 7si + ^7«i + 1 

where 5 S and S = {Si\i E Mc} are normalized transmission power of source and the set of transmission powers of 
relays, respectively. Here, S normalizes the transmission power to the case where no power adaptation is employed. 



For the sake of tractability of theoretical performance analysis in some cases, we may use an upper bound on 7$ as 

7»,u = min(^7 si , -^7«i)- (6) 

This yields an upper bound on the equivalent SNR at the destination as 7 eq u = ^-"fsd + J2ieM c ^> u ' ^his i s a 
meaningful bound frequently used in the literature [5|-|6| and facilitates the analysis in many cases. 

We use a set of N transmission modes (TM) each corresponding to a combination of modulation and coding. 
These TMs provide transmission rates of R\, R2, Rn bits per symbol and we assume R n > R n ~i- We can 
express the BER of TM n as an approximated function of received SNR, 7, through curve fitting by 

0.5e-f»T* n if 7 < 7 n 



BER n ( 7 ) 




(l + e c„( 7 -6»))' 

where p n , q n , a n , b n , c n and k n are the approximation fitting parameters. The value of 7^ is determined by the 
intersection of the two parts of the approximation. The expression for the second part is a modification of what is 
proposed in [ 1 1 1, and the first part is devised for better performance modeling in the low SNR regime. The inverse 
of BER„(7) describing received SNR as a function of BER is given by 

if P e > BER n ( 7 ^) 

(7) 

b n if P e < BER n ( 7 fJ. 

The objective is to maximize the spectral efficiency of the system while providing both security against eaves- 
dropping of untrusted relays and reliable communications for the destination. These two requirements are expressed 
as follows 

CI. BER, > BER^g, fori < i < N R 

C2. BERj < BER^t, (8) 

where BER; is the BER at the UR; node and BERj is the BER at the D node applying MRC. The average spectral 
efficiency is expressed as 

N R 

,, = ^f Pr(TM = n), (9) 

n=l 

where Pr(TM = n) = P n is the probability of selecting TM number n for transmission. The factor 1/2 multiplied 
by R n is due to the half-duplex cooperative transmission. The equation (O indicates that for calculation of average 
spectral efficiency of the system, we only need to derive an expression for P n . 



III. Link Adaptation and Untrusted Relay Assignment for Cooperative Communications 
The general optimization problem for LAURA with power adaptation is expressed as follows 

max max max R 

Mc T»G{l,2,...,iV} S S ,S 

s.t. CI. BERj > BER^for 1 < j < N R 
C2. BER(j < BERjJ^ 

C3. S s + S > ^ S tot, (10) 

iGMc 

where R is the transmission rate in a given frame in bits per symbol or equivalently the spectral efficiency. Note that 
the solution to this problem provides the instantaneous (per frame) power allocation scheme in the two transmission 
phases. A special case of this problem with Nq = 1 is considered in fT[, where a single relay is selected for 
cooperation out of the Nr available relays. In the following, the exact solution to the problem ( TTOb is presented 
and in the next section we explore some other possible solutions. 

A. Exact Solution 

The problem ( fTOb can be solved exactly and the resulting solution is referred to as LAURA1 in the sequel. The 
following proposition enables the proposed solution. 

Proposition 1: Identifying the TM with maximum R n and satisfying C2 in ( [Tol l, is equivalent to maximizing 

7eq- 

Proof: Consider 7* as the optimized 7 eq constrained to C2. Suppose that j' eq < 7* q . Now, if 
ri = argmax„{BER n (7^ q ) < BER? gt } and n* = argmax„{BER„(7* q ) < BER? gt }, then R n * > R n ,. ■ 
Using the proposition 1, for a given set A4c, the design problem ( TTOb may be expressed as follows 

max 7eq (11) 

s.t. CI. 5 S < 

C3. >S S + ^ ' Si < St t, 

ieMc 

where i — argmax^ 7^ and based on <j7j, = r„(BER[ gt ). Once this subproblem is solved, we can optimize the 
resulting solution over the choice of relay set, Aic- It can be easily shown that 7 eq is a concave function of (5 S , S). 



In other words, the next three conditions for 7^ hold E21 Appendix 1]: (d 2 ^ /dSf) < 0, (d 2 ^ /dS 2 ) < and 
(<9 2 7e ( q } /dSf)(d 2 ^ /dS 2 ) ~ [(d 2 j i ^/dS s dS l )} 2 < 0, for i e Mr, where 7^ denotes 7eq considering M c as the 
hh subset of A4r. Thus, the KKT condition gives the optimal solution to this problem. 

In order to simplify the solution, we first solve ( fTTT l without CI. The governing Lagrangian for this problem 
temporarily ignoring CI is 

C = 7eq + Ai ( 5 S + S i ~ 5 ™ ) ■ (1 V 

V iEMc / 



The optimal S s and Si should satisfy 



Then, 



that yields 



|J = and |^=0 foneMc. (13) 

dS s oSi 



dSi-Js-—^, ^2+^1-0. (14) 



S t = — 



f 7s* + if 7*d + l)' 



1 (S 2 ^ 2 , 5 S \ 1/2 S. 
v 



^3 7si7*d + ^2 7 si7id I - -<77s-. 



(15) 



where f = \/— Ai and denotes max(a:,0). Also, 

2 

7s»7id c2 1 TsiTid o 

2-, f a 7 . + a Td + i) 2 +Al ~°- (16) 

Substitution of ( fT5T l in C3 of (fTTT l and ( fToT ) gives a set of two equations. Numerically solving this set of equations 
yields the optimal Si for i E A4c and also S s . Now we check whether the solution satisfies CI in ( fTTT> . If the 
condition is satisfied, then the solution is the optimal one, else we set S s — S^p- in (I15t . Next, substitution of (fT~5T > 
in C3 of ( TTTb yields the optimum £. 

Using the presented solution to (fTTT l. Algorithm 1 formulates the LAURA 1 scheme or the exact solution to design 
problem (TlOb . 

B. Analytical Results for Nc = 1 

In case only one relay is to be selected for cooperation, i.e, Nc = 1, closed form solutions for transmission 
powers of source and relay can be obtained. This in turn allows for an analytical performance evaluation in this case. 



_ -Bi-y/ef-^ 



For Nq = 1 and using UR,; node for cooperation, we have Si = S (ot — <S S and ( TTST l yields S sni = 
where 



HPi 
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Algorithm 1 : Exact Solution for LAURA 

1) Select n = N. 

2) If n = then go to outage mode and exit. 

3) Build the set A4a comprised of all Nc member subsets of A^r and index its members as Aic,i, 
1 <= /1 9 Nr [ 1 

l fc \±, 6, (iv R -iVc)!JVc! J ' 

4) For i = 1 to {NR N £), m 

4-1) Set M c = Mcj 

4-2) Solve CEfl and obtain S^ and i G M c 
4-3) Calculate 7^. 

5) Select the set of cooperating relays by I* = argmax; 7^ and set 7 eq = 7^ \ 

6) If C2 in ( [Tol l is satisfied, set TM = n and exit; else n = n — 1 and go to step 2. 



=lfd(%i + 7sd) + 7s 2 i(7sd ~ lid) - 27 sd 7id7si 
0i = 2[7 S j(7 S d -7id +7id7sd<Stot) - 7 l 2 d 5 'tot(7s i +7sd)-7id7s< 



Pi = 7sd + 7 l 2 d(7si + 7sd)S' t 2 0t + (27 ld 7sd + 7id7si)<St c 



Then, the optimal transmission power of S node cooperating with UR S ; node is 



7s i 

The constraint C2 in ( TTOb in this setting is expressed by the following event 



A* ■ 7 (l) 



eq 



c 7sd 



. Sal—St 



-ltd + 1 



> r d 



J S S =S"' 



(17) 



(18) 



where T„ = r„(BER, gt ). For the presented transmission strategy in Algorithm 1, the event that TM number n 
satisfies both security and reliability requirements is denoted by A n and its probability is 

Pr(A„) = Pr ( U 4 J = 1 - Pr f Q Wj = 1 - E 7.,7 s5 |Pr Q W IT*, 7 S H | , 



1-E 



7s<i,7 s ; 



n[!- Pr (7e?>r?j7sd,7 s ;, 



(19) 



where E 7sd 7s , denotes the expectation with respect to independent variables 7 s d and 7^. The last equality in $1% 
results from the independence of A\ and A 3 n for i =^ j and given j s a and j^. The probability of selecting TM to is 

Pm = Pr ( (J A n j - Pr ( |J A„ J = E {I (A» m ) - I (A Wl ) } , (20) 

\n—m / \n— m+1 / 

where I(-) is the indicator function that is 1(E) = 1 if E is true and = if i£ is false, and 

*U7l?,7rf) = arg nun (7$ > F d „ | S s = S^O . (21) 

In fact, using UR; relay node with optimum power allocations described, nmileq , %{) denotes the minimum TM 
number from the set {to, N} for which the reliability and security constraints in given realizations of the fading 
channels are satisfied. Here, in order to facilitate the analysis we introduce an approximation whose effectiveness 
is shown in Section [VT] We consider 

M7^,7rf)wm. (22) 



Then, the probability of TM to is 



P ro = Pr(An)-Pr(>Wi)- (23) 



Appendix [A] presents the joint PDF of 7^ and 7 S j that is required to compute P m in d23l ) and hence the average 
spectral efficiency. The authors also presented the above analysis in (TJ. 

IV. LAURA: Constant Power Strategies 

The exact solution for ( fTOb is optimal and provides a benchmark for comparison to other possible simpler schemes. 
In the following, two suboptimal power allocation strategies are introduced. 

A. Constant Power Transmission (LAURA1-CP) 

A constant power solution to dTOb is obtained by setting = S and S s — S and removing step 4-2 in Algorithm 
1. In this case, the SNR of S-URi and 7 eq are required at the S node as CSI to enforce the security and reliability 
constraints, respectively. The equivalent SNR, j eq , is estimated at the destination and fed back to the S node, or it 
can be calculated by S node knowing CSI of S-D, S-UR^ and UR^-D. 

The performance of LAURA 1 -CP is evaluated in Section [VI] As elaborated below, for the case with one 

(i) 

cooperating relay (A*c = 1), the theoretical performance analysis is possible. Consider i* — argmaxi 7 eq obtained 
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from step 5 of Algorithm 1. The event that relay i* satisfies the reliability constraint corresponds to the event that 
there is at least one relay that can satisfy this constraint. For LAURA 1 -CP with single cooperating relay, we have 

N R N R 

A n : p| 7s, < r r „ n |J 7 « > It (24) 



(.=1 



i=l 



Then, the probability of selecting TM m is 



/ N \ /( JV R jv r -\ n r \ 

p£ p =p r \A m nf| Al = Pr U f| 7si < r m n |J 7 « > r d m n f) 7 « < r d m+1 

\ n=m+l / VI i=l i=l J i = l / 

/ r JVr -\ N R \ 

=Pr N H 7s* < n 7e ? < jn |J 7 w > r d J 

/JVr \ /JVr \ 

=Pr f) 7si < P r m n 7 « <r d „+i -Pr f] 7si < P r m H 7e ? < P d m . 



(25) 



Vi=l 



Since different relay channels are independent, given 7 s d, 7^ is independent of 7^ for i ^ j. Then 

JVr JVr 

pg= e 7 .{ n pr ^ r ™ n 7^ < r d m+ ii7sd) n pr (7- < r ™ n 7 e ? < r d m i 7 sd) }. 



(26) 



i=l i=l 

Using the upper bound for equivalent SNR in ©, a closed form expression for the TM probability is obtained as 



follows 



( JVr 

^„ p =Ejn|i-< 




e i« 



-e ->•<■ 



(27) 



Then, the average spectral efficiency is computed by substitution of (l27l i in ([9). The authors also reported the above 
analysis in (T|. 



S. Constant Power Relay Transmission (LAURA1-CPR) 

The LAURA 1-CPR scheme involves adaptive power transmission for the source and constant power relay 
transmission. In this case, we set Si = S for i £ Aic- The required CSI are S-URi and 7 eq . For 7 eq , we need to 
transmit from S node with the power S s , and wait for the estimation of 7 eq at the destination or it can be calculated 
by S node knowing CSI of S-D, S-UR^ and UR^-D. A modification of the Algorithm 1 yields the solution for 
constant power relay transmission, where, step 4-2 is replaced as follows 



Set Sf ) = S and S„W = min(S tot - N C S, S^-). 



7 S : 



(28) 
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V. LAURA: Relay Selection Strategies 

The solution presented for problem ( TTOb in Algorithm 1 provides high performance but involves an exhaustive 
search over all Nc member subsets of A4r. Specifically, it examines optimized jeq for the said relay subsets, which 
of course only a single one is finally used. In addition, the source requires jid for non-cooperating relays as CSI as 
well. Suboptimal solutions may have some advantages in practice, as their CSI requirements and complexity could 
be far less than the optimal one. In this section, we propose two efficient relay selection strategies for LAURA. 

A. Relay Selection Based on Source Relay CSI (LAURA2) 

An efficient solution with manageable CSI requirement may be constructed by taking a suboptimal relay selection 
approach that relies on 7^, i £ A4r and j eq resulting from Aic- Hence, there is no need to URrD CSI feedback 
(i G .Mr). To limit the complexity, we wish to avoid solving ( fTTT i to obtain 7 eq © for all subsets of VWr. To 
this end, the instantaneous source to relay SNR can be used as the relay selection criterion. Due to the security 
constraint, 7, is limited to a . As a result, the subset of relays that provide a high 7 eq may also be identified by 
selecting the relays with the highest 7s»'s. Indeed, as we shall see in Section [VlJ this suboptimal and yet efficient 
relay selection strategy does not significantly degrade the performance of the optimal solution. Algorithm 2, gives 
this solution that is labeled as LAURA2. 

B. Constant Power Relay Transmission with Relay Selection Based on Source Relay CSI (LAURA2-CPR) 

The relay selection criterion according to SNR of source to relays can also be applied to the LAURA1 scheme with 
constant power relays (Section [IV-Bl i. This will further reduce the computational complexity and CSI requirements. 
A modification of the Algorithm 2 yields the solution with constant power relay transmission, where the step 4 is 
replaced with the following 

Set S Z = S and S s = min(S m - N C S, S^). 
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Algorithm 2: LAURA with Modified Relay Selection 

1) Select n = N. 

2) If n = then go to outage mode and exit. 

3) Sort the set of available relays .Mr according to their7 S i in descending order and select A4c as the first Ac 
in the set. 

4) Solve (flU and obtain S s and S» for i e Mc- 

5) Calculate 7 eq . 

6) If C2 in (TTOb is satisfied, set TM = n and exit; else n = n — 1 and go to step 2. 

C. Constant Power Relay Transmission with Relay Selection Based on Source Relay Channel Statistics (LAURA3- 
CPR) 

The instantaneous relay selection criterion according to SNR of source to relays involves high speed switching 
of relays per frame and hence a rather sizable network control overhead. A relay selection criterion according to 
the statistics of source to relay channels can help mitigate this problem. This will further reduce the computational 
complexity and CSI requirements since the selected relays are fixed as long as the average SNRs of source to relays 
remain unchanged. A modification of the Algorithm 2 yields the solution for constant power relay transmission in 
which the step 3 is replaced as follows 

Sort the set of available relays .Mr according to their 7^ in descending order and select Aic as the first Ac in the 
set. 

The selection of relays according to the average channel conditions also makes the theoretical performance 
analysis of the system possible. In the following, we present a performance analysis of LAURA3-CPR. The event 
A n in this case using the upper bound for equivalent SNR is expressed by 

Pr(A„) = Pr ( 7eq , u >T d n \S s = § s ,„( 7s i)) , (29) 

where <S SjTl (7 s ^) = min(S' tot — NqS, S^t). In order to calculate the probability of event A n , we first derive the 
moment generating function (MGF) of 7 eq , u given 7^ and S s — S s ^ n (j s i)- Since |f7sd and {7i, u } are all independent 



of each other given 7 j, the desired MGF knowing 5 S = S , s ,n(7 s i) * s gi ven by 



M^ uhsl (s\x) =M^ h ,(s\x) Af 7<i ,| 7i j(«|aj), (30) 

ieMc 

where Ms^^^a;) and M 7ill | 7s ,(s|a;) are respectively the MGF of 4f7sd and 7; jU given 7 s j and S s = ^ s ,n(7 s -). 
Using the definition of the MGF as Mx(s) = E (e _sX ), it can be easily shown that 

M Ss , (six) = —~ . (31) 

In order to calculate M 7j j 7 ,(s|a;) we first calculate cumulative distribution function (CDF) of as follows 

F 7iu | 7sJ (z) = Pr( 7l , u < z\ ls , =x) = l- Pr(^ 7sl > z\ lA = x)Pr( 7ld > z), (32) 
where Pr(7id > z) = e^J. Using the conditional PDF of f-y si \-y,(y\x) according to Appendix lAl we have 



S 

Pr(^7si > z\%i = x) 



S s ,nM/S 



Bi(x) 



gj(g) 1 r -3h u < T , A 1 Di W fi( v r ) 
C(z) 7s, y ) + CW %_I) 



5; 



Ss,n(x) 



(33) 



C(x) V / C(x) 

where U{ ) and <5(-) are unit step and unit impulse functions, respectively, and Bi(x), C(x) and Di(x) are defined 
in Appendix lAl Then, 



M 7 . u | 7r (s|a;) = s / F lt u | 7 ,(z|a;)e sz dz 
Jo 

. Bi(x) s 



C(x) 

fB^x) =2. Di(x) 



1 - e 



S.,n(s0 

-a— — 



Ss,„(*Hsi 7id 



lew 



C(a:) / ^ + s 

V y 7 7»d 



1 - e" 



.gjtnW ( 1 



Finally, 



Pr(A n ) = Pr ( 7eq > r*) = 1 - E x {Pr ( 7eq < T d „) | 7§J = x} = 1 - E J £" 



-1 f^KiC 5 !^) 



(34) 



(35) 



7«,=r d „ 



where £ _1 denotes inverse Laplace transform with respect to 7 equ that is simply computed through symbolic 
evaluation with MATLAB for every AA R and Mq- Then, the expectation with respect to 7 s j is computed through 
numerical integration. 

The probability of TM m for LAURA3-CPR is expressed as 



P. 



LAURA3-CPR 



Pr ( |J A n j - Pr ( |J A n \ - E {I (A hm ) - I (A Am+1 ) } 



(36) 
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where 



n m (7eq, u , 7 s j) = arg min 



n£{m,...,N} 



(7 eq ,u>r*|5 s = 5 s , n ( 7sl )) 



(37) 



In order to facilitate the analysis we consider an approximation like the one in d22i > as 



«m(7eq,u,7si) 



m. 



(38) 



And finally, the probability of TM m is 



P 1 



,LAURA3-CPR 



Pr(A ro )-Pr (A m+1 ). 



(39) 



m 



Calculating the probability of TM based on d39l and (1351 . enables the computation of average spectral efficiency 
according to ©. 

Remark 1: Since for iVc = A/r no relay selection is employed, the LAURA performance is independent of the 
relay selection strategy. Hence, the presented performance evaluation of LAURA3-CPR in d39l (and d35l l) also 
applies to LAURA1-CPR and LAURA2-CPR schemes. 



In this section, the performance of the proposed LAURA 1, LAURA2 and LAURA3 schemes of different relay 
selection strategies in conjunction with different power control mechanisms of adaptive power, CP and CPR 
are evaluated. Both analytical and numerical results are presented and the effects of different parameters on the 
performance are investigated. 

A. Experiment Setup 

Figure [T]illustrates the topology of the network under consideration with A/r available relays. We consider average 
SNR of each channel proportional to due to path loss, where I denotes the distance between the two parties (a 
is set to 4 in this paper). Without loss of generality we consider the distance between S and D nodes normalized 
to 1. Relays are all located on a line perpendicular to the line connecting S and D nodes each distanced 0.1 apart 
as depicted in Figure Q] As shown in Fig. Q] URi node is fixed at the distance 4i on the line connecting S and 
D nodes. In all figures, BER[ gt and BERf gt are set to 0.1 and 10 -6 , respectively. The sum power constraint S tot is 
set to (N-r + l)S in all cases. 

The TMs used with the proposed LAURA schemes could be in fact any set of possible channel coding and 
modulation pairs. The suggestion is to use sharp channel codes to obtain sharp BER curves for TMs, which provide 
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Fig. 2. BER curves for six transmission modes and their corresponding rates 

TABLE I 

Transmission Modes for AMC Scheme and Their Corresponding Fitting Parameters 



Mode 


1 


2 


3 


4 


5 


6 


modulation 


BPSK 


4-QAM 


8-QAM 


8-QAM 


16-QAM 


32-QAM 


coding rate 


1/2 


1/2 


1/2 


2/3 


3/4 


4/5 


rate:i? n 


1/2 


1 


3/2 


2 


3 


4 


Pn 


2.97 


1.17 


0.80 


0.65 


0.46 


0.38 


q-n 


1.05 


0.68 


0.54 


0.59 


0.67 


0.58 


On 


1.55 


0.13 


0.14 


0.12 


0.06 


0.07 


bn 


0.62 


1.43 


3.48 


6.19 


10.37 


24.09 


c„ 


41.9 


34.90 


20.03 


3.52 


5.13 


2.21 


kn 


16.81 


104.70 


48.83 


75.99 


6.77 


6.50 



a low security gap and an acceptable reliability performance (for the design parameter ranges of interest). The set 
of LDPC codes in DVB-S2 standard lfl3l offer this characteristic. These channel codes in conjunction with different 
modulation schemes yield the BER curves depicted in Fig. [2] The TMs specifications and fitting parameters are 
presented in Table U 

B. Numerical Results 

Figure[3]depicts the end to end average spectral efficiency of LAURA1, LAURA1-CPR and LAURA1-CP schemes 
illustrating the effect of type of power control. The horizontal axis indicates 7 s d. It is observed that for LAURA1-CP, 
increasing average SNR of S-D link beyond a specific value leads to a reduced average spectral efficiency. This is 
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due to the limitation in number of TMs. If there is a TM that could still satisfy the security constraint in high SNRs, 
then the declination occurs later at a higher average SNR of S-D channel. Comparing the performance of LAURA1 
and LAURA 1 -CP highlights the advantage of source power adaptation that significantly reduces the outage events 
for enforcing the security constraint. 

The performance comparison of LAURA 1 and LAURA 1-CPR suggests that for low to medium SNR regimes 
relay power adaptation improves the performance of the system. This is more significant with larger Ac's. One sees 
from the performance of LAURA 1 and LAURA 1 -CP that the larger the number of cooperating relays, the better 
the performance we achieve through cooperation. However, there is no significant advantage in using more than 
four relays for Ar = 5. The same statement also holds true for LAURA1-CPR up to Ac = 4. This observation is in 
contrast to the understanding in traditional (non-secure) amplify and forward relaying that the best relay selection 
achieves almost all of the spectral efficiency performance gain fl4l . Comparing LAURA 1-CPR for Ac = 4 and 
5 in low SNR regimes reveals that with constant power relays and sum power constraint, increasing the number 
of cooperating relays beyond a limit may decrease the performance. The reason is that in this case, the available 
power to the S-node (and hence the SNRs of source to relays) becomes limited. 

Figure|4]illustrates the effect of relay selection strategies on the average spectral efficiency performance of LAURA 
schemes. The main observation is that the performance of relay selection according to source-relay channel SNRs 
coincides with that of optimal relay selection. This is evident both by performance comparison of LAURA 1 with 
LAURA2 and LAURA1-CPR with LAURA2-CPR. It is also observed that the performance of LAURA3-CPR with 
relay selection according to average source-channel SNRs approaches that of LAURA2-CPR utilizing instantaneous 
channel SNRs only when the number of relays increases to Ac = 5. Nevertheless, the relay selection according to 
average channel statistics offers softer relay switching. 

Figure [5] depicts and compares the simulation and approximate analytical results for certain LAURA schemes. As 
discussed in Remark 1, for Ac = Ar, since no relay selection is employed, the theoretical calculation of spectral 
efficiency proposed for IV-CI applies to all the schemes with constant power relay transmission. The results verify 
the accuracy of the proposed analytical performance evaluations. 

Figure [6] demonstrates the effect of relays positions on the performance of the proposed LAURA scheme. It is 
evident that when the relays are positioned closer to the D node, LAURA provides a better average spectral efficiency 
performance. In fact, in this setting the security constraint is satisfied more easily. Our experiments reveal that in an 



amplify and forward cooperative communication system, in the SNR ranges of interest for cooperation (low S-D SNR 
regimes), when a sufficiently large number of available relays are utilized, provisioning the security constraint in 
LAURA only imposes a negligible spectral efficiency performance penalty. In fact, for larger number of cooperating 
relays, the performance penalty is small even for high SNR regime, e.g., less than 7% for Nq = Nr = 5. 

VII. Conclusions 

In this paper, a link adaptation and untrusted relay assignment framework for cooperative communications with 
physical layer security is proposed. The design problem is set up for highly spectrally efficient communications with 
reliability for the destination and security in the presence of untrusted relays. The security constraint is imposed 
by ensuring that the relays cannot decode useful information from the signal they relay. The optimal solution to 
the design optimization problem is presented while strategies for mitigating practical challenges are also proposed. 
This involves several relay selection strategies and power adaptation solutions. Performance of these approaches is 
analyzed theoretically, in certain cases, and rigorously through simulations. The effect of different design strategies 
and parameters, including relay selection, power control at source and/or relays, relays positions, and number of 
cooperating relays are investigated. 

The future research of interest in this direction includes considering the (non-cooperative) behavior of relays 
without a service level trust. One example is when the relays use MRC in contrast to the agreed relaying protocol, 
where they may be treated as eavesdroppers. An alternative example is when the relays may feed back untrue CSI. 
Another orientation of research interest is to design particular channel codes with low security gap and reliability 
performance within this framework and in line with the works reported in f9j- 

Appendix A 

In this section we derive the joint and conditional PDF of 73 and 7^. The joint cumulative distribution function 
of these two variables is given by 

N R 

F^siix^y) =Pr(max7 S j < x,^ % < y) = JJ ( 1 - e~w) - e~^U(y - x) + (l ~ e~^ijU{x - y) 

3& 
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Then, the joint probability density function of 7 s j and 7^ is 





dxdy 
Nr 1 

1=1 &i 



U{x - y) 



+ 




i** S(y — x) 



=Bi(s) — e *«U(x-y) + Di(x)8{y - x), 



.7si 



(40) 



and the PDF of 7^ is simply 



N R 1 



(41) 



The PDF of 7 si given 7 s j is then 



hi\iM x ) = 



hiTM ( x >y) = Bj(x) [J_ e -^- 

f %J (x) C{x) [7/ 



U(x-y) + 



Di(x) 



C(x) 



8{y-x). 



(42) 
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Fig. 4. Average spectral efficiency of different relay selection schemes for Nr = 5 and £ s i = 0.9. Solid, dashed and dotted lines are for 
LAURA 1, LAURA2-CPR and LAURA3-CPR, respectively. No line is used to connect the markers associated with LAURA2, but they closely 
follow those of LAURA 1. Note that the performance of LAURA2-CPR and LAURA3-CPR for N c = 5 coincide. 
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Fig. 5. Comparison of analytical and simulation results for average spectral efficiency of different LAURA schemes for A/r = 5, 4l = 0.9. 
Solid and dashed lines depict the simulation and analytical results, respectively. 
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Fig. 6. Average spectral efficiency of LAURA 1 scheme for different relay placements. 



